Inequivalence of pure state ensembles for open quantum systems: 
the preferred ensembles are those that are physically realizable 
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An open quantum system in steady state p ss can be represented by a weighted ensemble of pure 
states p ss = pfc|^>fc)(^>fc| in infinitely many ways. A physically realizable (PR) ensemble is onc 
for which some continuous measurement of the environment will collapse the system into a pure 
state \tj}(t)}, stochastically evolving such that the proportion of time for which \t[>(t)) = \ipk) equals 
pfe. Some, but not all, ensembles are PR. This constitutes the preferred ensemble fact. We present 
the necessary and suflicient conditions for a given ensemble to be PR, and illustrate the method by 
showing that the coherent state ensemble is not PR for an atom làser. 
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Introduction. It is well known that there are infinitely 
many ways to write a given mixed state as a mixture 
of (in general non-orthogonaí) pure states. That is, 
there are infinitely many ensembles, {(n^, pfe)}&, consist- 
ing of ordered pairs of pure states (rank-one projectors 
íïfe = |V ; fc)(V , fe|) an( A weights pfc, that represent a given 
state matrix (density operator) p according to 



^2 P^fc. 
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Any such ensemble representation suggests an ignorance 
interpretation of the state matrix p. That is, one would 
claim that the system "really" is in one of the pure states 
ITfc, but that one happens to be ignorant of which state 
ITfc (i.e. which k) pertains. The weight p^ is interpreted 
as the probability that the system has state ílfe. 

This interpretation can only be maintained for so- 
called proper mixtures, i.e. those for which the system 
is not entangled with its environment [jïj. However, an 
improper mixture may be turned into a proper mixture 
simply by appropriately measuring the environment, and 
ignoring the result. The partition ensemble fallacy Q or 
preferred ensemble fallacy is that one may use a partic- 
ular ensemble to draw inferences about an experiment, 
even if those inferences depend upon the choice of ensem- 
ble. Here we do not deny that this is indeed a fallacy. 

Mixed states arise naturally in the description of open 
quantum systems, in fields as diverse as chemical physics, 
quantum òptics, micromechanics and nano-electronics. It 
is often possible to approximate the evolution of such sys- 
tems by a Markovian master equation p = Cp, where C 
is the Lindbladian j|. Under these conditions, an ex- 
perimenter may perform continual measurements on the 
environment with which the system interacts without af- 
fecting the master equation evolution (HJJ]. Hence the 
system state can be considered a proper mixture. 

Let us restrict the discussion to a stationary mixed 
state p ss , assumed unique: 



0. 
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Also, let us consider only stationary ensembles for p ss . 
Clearly, once the system has reached steady-state then 
such a stationary ensemble will represent the system for 
all times t. Now, if the ignorance interpretation were to 
hold for such an ensemble then it should be possible, in 
principle, for the experimenter to know which state Tlk 
pertains to the system at any particular time t. The per- 
tinent index k would change stochastically such that the 
proportion of time the system has state Tík is pfc. 

In this letter we show that for some ensembles this 
ignorance interpretation fails. That is, there is no way 
that an experimenter can know at all (long) times that 
the system has some ensemble state tlk- We say that 
such ensembles are not physically realizable (PR). How- 
ever, there are other stationary ensembles that are PR. 

The existence, for a given system, of two non-empty 
classes of stationary ensembles, those that are PR and 
those that are not, constitutes a preferred ensemble fact 
(PE-fact). Note that this fact identifies a preferred class 
of ensembles (the PR ones), not a unique preferred en- 
semble. Also note that this fact does not contradict the 
argument against preferred ensembles in Ref. Q , because 
it is a fact about ensembles representing a stationary 
state for all times. As we will see, the PE-fact has sur- 
prising implications for some open quantum systems. 

This letter is organized as follows. First we discuss the 
Hughston-Josza-Wootters theorem || and why it does 
not contradict the PE-fact. Then we give the necessary 
and sufficient conditions for an ensemble to be PR for a 
given system. For linear systems with uniform Gaussian 
ensembles we show these conditions reduce to a simple 
inequality. Next we establish the PE-fact for a particu- 
lar system of interest: a model atom làser. Finally we 
discuss the implications of our results. 

The HJW Theorem. We wish to consider a system with 
state matrix p which is mixed solely due to its entangle- 
ment with a second system, the environment. That is, 
there is always a pure state |^) in a larger Hilbert space 
of system plus environment such that 
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Certainly it is always possible formally to find such a 
pure state [^J. Also, if the environment is initially not 
pure then it can be then measured in its diagonal ba- 
sis so as to make it conditionally pure without affecting 
its subsequent interaction with the system. Therefore we 
can assume the existence of |\E r ) without loss of generality. 

The Hughston- Jozsa-Wootters (HJW) theorem || says 
the following. For any ensemble {(Tlk, pk)}k that repre- 
sents p, it is possible to measure the environment such 
that the system state is collapsed into one of the pure 
states IT/j with probability pk- That is, there exists a 
probability-operator measure {Fk}k acting on the envi- 
ronment Hilbert space such that 



pktík = Tr c 



sys 
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The HJW theorem gives rigorous meaning to the igno- 
rance interpretation of any particular ensemble. It says 
that there will be a way to perform a measurement on 
the environment, without on average disturbing the sys- 
tem state, to obtain exactly the information as to which 
state the system is "really" in. Of course the fact that 
one can do this for any ensemble means that no ensemble 
can be fundamentally preferred over any other one, as a 
representation of p at some particular time t. 

The HJW theorem does not contradict the preferred 
ensemble fact introduced above. This is because the PE- 
fact refers to the much stronger notion of representing the 
state of the system obeying a master equation, at steady 
state, by a stationary ensemble that applies at all times. 
This extra condition means that the HJW theorem does 
not apply, and not all ensembles are PR through mea- 
suring the environment of the system. We will now give 
the conditions for an ensemble to be PR. 

Establishing the Preferred Ensemble Fact. Consider a 
system obeying the master equation p = Cp in the long- 
time limit so that it is described by the (assumed unique) 
stationary state p ss . Now say we wish to give an igno- 
rance interpretation at all long times to a particular sta- 
tionary ensemble {(ÍI^, pk)}k satisfying 



^Pk% 
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At a particular time t this is always possible by the HJW 
theorem. That is, we can measure the environment to 
find out the pertinent state l·lfc, and, on average, the sys- 
tem remains in state p ss . This may involve measuring 
parts of the environment that interacted with the sys- 
tem an arbitrarily long time ago, but there is nothing 
physically impossible in doing this. 

Now consider the future evolution of the system state 
following this measurement. At time t + r, it will have 
evolved to pk(t + r) = exp(£r)ÍTfc. This is a mixed state 
because the system has now become re-entangled with 



its environment. The system state can be repurified by 
making another measurement on its environment. How- 
ever, if the same ensemble is to remain as our represen- 
tation of the system state then the pure system states 
obtained as a result of this measurement at time t + r 
must be contained in {ÍI/j}?;. By the HJW theorem, this 
will possible if and only if pk (t + r) can be represented 
as a mixture of these states. That is, for all k there must 
exist a probability measure {u^fefr)} such that 



exp(£T)n fc = Wj k (T)ílj. 
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If Wjk(T) exists then it is the probability that the mea- 
surement at time t + r yields the state ÍL; . 

Equation (jfy is necessary but not sufficient for the en- 
semble {(£[?, pj)}j to be PR. We also require that the 
weights be stationary. That is, that for all j and all r, 



Pj 



k 



PkWjk(r). 
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Multiplying both sides of Eq. (|g) by pk, and summing 
over k, then using Eq. (0) and Eq. (g) gives e Cr p ss = p ss , 
as required from the definition of p ss . 

If the two conditions (0) and (0) are satisfied for some 
time t then they will be satisfied for any múltiple of r. 
Therefore it is sufficient to establish that they are sat- 
isfied for t = dt. (Here dt is infinitesimal with respect 
to all relevant system time scales but strictly must be 
long compared to the environment correlation time; the 
master equation is not vàlid on any shorter time scale.) 
For a bounded Lindbladian superoperator C and distinct 
states {l·lfcjfc, if Wjk(dt) exists then it is given by 



w jk {dt) = ó kj - 5 kj ^2 Ikidt + -fkjdt 



(8) 



with "fkj finit e. The coefficient 7^ may be interpreted 
as the rate for the system to jump from state Tlk to 
state Tlj. That is, the quantum master equation is re- 
placed by a classical master equation || for the occu- 
pation probabilities {pk}k f° r the states {l·lfcjfc where 
Pj(t + dt) = J2kPk(t)wjk{dt), and so from Eq. (0) 



Pk 
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ikj + ^ik 3 Pj- 
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The condition (Q) is then equivalent to the condition that 
Pk = pk be the stationary solution of Eq. (^|) . 

Any stationary scheme of continual rank-one measure- 
ments of the environment will produce stochastic dy- 
namics of this sort in the steady state j?]]. The igno- 
rance interpretation of the ensemble thus produced is 
then rigorously justified. The PE-fact follows if, for the 
given system, there exists at least one ensemble satisfy- 
ing Eq. (||) for which there does not exist a probability 
measure {wjk(r)}j satisfying conditions (g) and (^). 
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The replacement of a quantum master equation by a 
classical master equation where the classical index k is as- 
sociated with a quantum state fl^ is often used to provide 
an intuitive picture of irreversible quantum dynamics. A 
canonical example is Einstein's description Q of an atom 
driven by a thermal field in terms of jumps from ground 
to excited states (absorption) and from excited states to 
ground (emission). Of course Einstein did not know the 
more general description in terms of the quantum master 
equation, but it is easy to verify that the the ensemble 
consisting of atòmic energy eigenstates is a PR ensemble. 
Specifically, to realize this ensemble one must continually 
measure the environment in the photon number basis. 

Linear Dynamics and Uniform Gaussian Ensembles. 
The description given above applies most naturally to 
ensembles with a discrete set of states {ÍIfc}fe. In many 
cases we wish to consider a continuum of states. In these 
cases it is convenient to take a limit in which the jump 
process described by rates 7^ is replaced by a diffusion 
process. We restrict ourselves to systems with linear dy- 
namics, and uniform Gaussian ensembles. These terms 
(defined below) only apply to quantum systems whose 
state p can be represented by a Wigner function || on 
2rí-dimensional Euclidcan phase space: 



W{z) = Tr 



p s íï / exp 

m=l ^ 



(10) 



Here S denotes ordering symmetric in the operators 
0*1, • • • z 2n ) = (ii,pi,x 2 ,p 2 , • • • ,x„,p n ), 



(11) 



where the co-ordinates {x n } n and momenta {p n } n each 
form a mutually commuting set of operators with the 
reals as eigenvalues, but (with h = 2) [x ni p n >] — 2iS n ,n'- 
Such a system has linear dynamics if and only if (for a 
suitable choice of co-ordinates and momenta) its Wigner 
function obeys an Ornstein-Uhlenbeck equation |6| 



W(z) = (V 



DV) W(z). 



(12) 



Here V TO = d/dz m , K is a constant matrix, and D is a 
constant matrix that is symmetric and positive semidef- 
inite. For simplicity we will assume that the eigenvalues 
of K have positive real parts, so that the system has a 
stationary state 



W ss (z) = G(z;0,V ss ), 
where V ss is defined by Q 

KV SS + V SS K T = D. 



(13) 



(14) 



Here we are using the notation that G(z; fl, V) is a multi- 
variate Gaussian in z parameterized by the mean vector 
fl and the covariance matrix V m . 



A uniform Gaussian ensemble of pure states 
{(ílfc, pk)}k comprises states í\ k that have Wigner func- 
tions Wk{z) that are Gaussians 



W k (z) = W(z;p k )=G(z;fl k ,V) 



(15) 



with mean determined by k but variance independent of 
k. That is, the states all have the same "shape" , but have 
different displacements in 2n-dimensional phase space. 
Since the ensemble represents p ss we have (converting 
from the discrete k to the continuum variable fl) 



W ss (?) 



dpi 



dfj, 2n p(p)W(z;fl). 



(16) 



Since W ss (z) and W(z;fl) are Gaussians, p(fl) is also: 

p(/J)=G(M;0,V 8S -V). (17) 

Now consider the conditions for the uniform Gaussian 
ensemble to be PR, starting with Eq. (^|) . From Eq. ( |Ï2] ) , 
if the system begins with Wigner function W(z; /2) then 
after an infinitesimal time dt its Wigner function will be 

G(z; fi - Kpdt, V + T>dt- KVdt - VK T dt). (18) 

This can be writtcn as a mixturo of the ensemble states 



dp'w(/j, ; fl; dt)W(z; jl) 



(19) 



if and only if the transition probability density is 
w{p'-p-dt) = G(p';p-Kpdt,T)dt - KVdt - VK T dt). 

(20) 



Since the ensemble is specified completely by V, Eq. 
turns into the condition on V that 



B 



v 



D KV VK >0. 



(21) 



An equivalent condition has been considered by Diósi and 
Kiefer ]9| in a related context, but they did not make the 
connection with physical realizability and measurements. 

Equation ( ^Ò] ) implies that the quantum master equa- 
tion in the Wigner representation (|ï^) is equivalent to a 
stochastic process for the mean displacements fl of the 
uniform Gaussian states with Wigner functions W(z; jl). 
The probability distribution for fl is governed by 



p(p) = ( V T K/I+ ÍV T B V V p(p), 



(22) 



where here V TO = d/dp m . This is analogous to Eq. (g). 

To establish that condition (Q) is also satisfied we thus 
have to show that the stationary solution of Eq. ( p2| ) is 
p(fl) = p(p). From Eq. © we get @ 



(23) 



where U is defined by KU + UK T = By. Now using 
Eqs. @ and ^) we obtain 
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U = V ss - V. (24) 

From Eq. (|l7j ) it follows that p ss {p) = as desired. 

Example: The Atom Làser. We now apply the above 
formalism to the problem of an atom làser. We con- 
sider only a CW atom làser (wl·iich has not yet been re- 
alized), which would consist of a continuously damped 
and replcnishcd Bosc-Einstein condensate. The damping 
at rate k would produce a continuous beam of coherent 
àtoms [Ï0| ]. The simplest quantum model for such a de- 
vice uses a single mode description of the condensate, 
with annihilation operator à = (x + ip)/2. Assuming a 
Poissonian pump with rate np, the stationary state of 
the làser is JÏÏj,[n| 

Pss = f = f>^^>)<n|, (25) 

Here \^fpe 1 ^) is a coherent state (eigenstate of à) |ÏI| . 

Equation ( |25| ) shows two of the infinitude of diffcr- 
ent ensemble representations of p ss . The number state 
ensemble {|?T-)} n is clearly PR by making continual mca- 
surements of the environment in the atom number basis. 
On the other hand, to determine the status of the coher- 
ent state ensemble { | ^/pé 1 ^) } , , it is necessary to examine 
the dynamics of the atom làser. 

The simplest reasonable model for an atom làser comes 
from taking the Standard ideal optical làser model [ p~0|JÏ~2| 
and adding the interaction between condensate àtoms 
p3[ , governed by the Hamiltonian (Tik/ '4p)xà^ 2 à 2 , with 

nm J 

Here is the wavefunction for the condensate mode, 
and a s is the s-wave scattering length. This nonlinear 
Hamiltonian is difficult to treat exactly, so we linearize 
the Lindbladian evolution about a mean ficld (à) = ^JJi 
(as appropriate for considering the physical realizability 
of the coherent state). We have shown elsewherc jÏ3| that 
this results in linear quantum dynamics for the Wigncr 
function W(x,p) as in Eq. ( |Ï2| ) with 

*"(;s) iD "(s»!»)- (-) 

Here x and p represent amplitude and phase fluctuations 
respectively, and v > is the excess phase diffusion. 

To check if the coherent state ensemble is PR we sim- 
ply need to check the condition (^ï|). We find 

v-(j;)^Bv-.(4 a? „). m 

The matrix By is positive semidefinite only for x = 0. 
That is, the coherent state ensemble is not PR for any 
finite atòmic interaction energy. 



Discussion. In this letter we have introduced the nec- 
essary and sufRcient conditions for an ensemble of pure 
states to be physically realizable (PR) as a stationary de- 
scription of a Markovian system in the long-time limit. 
Here physically realizable means that, by measuring the 
environment of the system, an experimenter could col- 
lapse the system state into a stochastically evolving pure 
state, such that the proportion of time it is in a partic- 
ular pure state is equal to the weight of that pure state 
in the ensemble. For uniform Gaussian ensembles for 
systems with linear dynamics we derive a simple inequal- 
ity to distinguish PR ensembles from non-PR ensembles. 
The existence in general of these non-empty classes of 
ensembles constitutes the preferred ensemble fact. 

To illustrate this fact, we have shown that in a simple 
model for an atom làser, the number state ensemble is 
preferred over the coherent state ensemble because only 
the former is PR. This is due to the presence of atòmic 
interactions, described by a Hamiltonian proportional to 
the a) 2 a 2 . We proved this using a linearized analysis 
which included the loss and gain processes. However, 
the result is perhaps not surprising since it is wcll known 
in quantum òptics that this Hamiltonian will turn a co- 
herent state into a squeezed state which cannot be 
described by mixture of coherent states. 

This result is interesting because the ignorance inter- 
pretation of the coherent state representation ( "the làser 
is really in a coherent state, I just haven't been bothered 
to find out which one it is") is very commonly adopted 
in quantum òptics (see Ref. for a discussion). Al- 
though it is tcnablc in principle for most optical làsers 
where x <C 1, it is not tenable for atom làsers where we 
expect x ^ 1 [0- So although atom làsers may give a 
coherent output, the condensate cannot be meaningfully 
considered to be in a coherent state. The relation be- 
tween output coherence, and the physical realizability of 
states of the condensate, will be explored elsewhere. 

HMW gratefully acknowledges conversations with T. 
Rudolph and R. Spekkens. 
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